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Abstract 

Recent experimental results from RHIC and LHC on hard photon emission 
rates in heavy-ion collisions indicate a large azimuthal asymmetry of photon emis- 
sion rate parameterized by the elliptic flow V2- Motivated by a recent proposal that 
the early magnetic field created by two colliding heavy-ions may be responsible for 
this large azimuthal asymmetry of photon emission rate, we compute the azimuthal 
dependence of the photon emission rate from a strongly coupled finite temperature 
plasma with magnetic field in the framework of gauge/gravity correspondence. We 
also propose and compute a new observable, "in/out-plane polarization asymmetry", 
constructed from the polarization dependence of the photon emission rates. We ob- 
serve that both the azimuthal and polarization asymmetry of photon emissions are 
strongly affected by the triangle anomaly (chiral anomaly) for low frequency regime 
below 1 GeV. 



*e-niail: hyee@uic.edu 



1 Introduction 



The quark-gluon plasma created in relativistic heavy-ion colhsions is an interesting new 
state of matter where one can test various theoretical ideas of QCD. The challenge is that 
this plasma lives for very short period of 10 fm, and there are many indications that it is 
strongly coupled, making perturbative computations very difficult. 

Their rapid evolution in time after about 1 fm is believed to be described fairly well 
by viscous hydrodynamics, but the earlier dynamics still remains as an important open 
problem. Once hydrodynamics sets in, the data needed to describe the system reduces 
dramatically to a few transport coefficients, and this is both a power and limitation of 
hydrodynamics. As the system is very likely strongly coupled, there has been much 
research on using gauge/gravity duality or AdS/CFT correspondence to describe the 
hydrodynamic regime of strongly coupled plasma, essentially computing the transport 
coefficients non-perturbatively [1]. Using gauge/gravity duality for early-time dynamics 
has also been pursued while it requires some amount of numerical analysis to solve relevant 
non-linear Einstein equations [21 El IH El E] . 

Some experimental probes see the created quark-gluon plasma in more refined way 
than the hydrodynamics description. A good example is the photon emission spectra 
from the plasma, especially at high frequency (and momentum) whose scale is above 
the temperature scale of the plasma. In linear response framework which is suitable for 
electromagnetic probes such as photon emission, the emission spectra mirror the Green's 
functions of vector current that the electromagnetism couples to. These Green's functions 
at high frequency-momentum are properties of the plasma beyond hydrodynamics, and 
they do depend on the microscopic details of the theory. 

This break-down of universality that has been acclaimed for hydrodynamics is an im- 
portant aspect one has to bear in mind in applying AdS/CFT correspondence to photon 
emissions, because for example N=4 Super Yang-Mills is very different from QCD micro- 
scopically. There is however another sense of universality that may originate just from the 
strongly coupled nature, and this still remains as a useful argument for applying AdS/CFT 
correspondence even for probes beyond hydrodynamics such as photon emission. 

Once photons are emitted from the plasma, they rarely interact with the plasma again 
because of the weakness of electromagnetic coupling. This aspect makes the photons ideal 
probes to the early dynamics of heavy-ion collisions, and there have been important exper- 
imental results from PHENIX [7] and ALICE [8J on photon emissions. At present, the 
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amount of photons observed experimentally tends to exceed the theoretical predictions, 
but theoretical computations are being refined to be closer to the experiments P, ITO] . 
What is challenging is that the observed photons (both direct and indirect) are from the 
entire history of the heavy-ion collisions, although different regions of frequency should 
be sensitive to different parts of the history. The low frequency regime will naturally be 
sensitive to somewhat later stage of plasma evolution with lower temperatures as well as 
hadronic regime after-burner phase. The higher frequency spectrum is expected to get 
dominant contributions from early stages, but the dynamics before 1 fm is still not very 
well understood. 

One striking recent result from PHENIX [7] and ALICE [8J is the large elliptic flow 
of direct photons at relatively large transverse momentum > 1 GeV. At = 2 GeV, it 
is as large as V2 = 0.25. Contrary to hadronic elliptic flow which builds up its magnitude 
through collective hydrodynamic evolution, high frequency photons coming from early 
stage are not expected to be sensitive to this late time collective evolution. (The low 
frequency spectra do get some effects from the Doppler shifts.) Given this, it is natural 
to seek a more direct source of azimuthal asymmetry relevant for the photon emission at 
the early stage of the heavy-ion collision. 

An interesting candidate is the magnetic fleld created by two fast moving colliding 
heavy-charged nuclei pXj- Its magnitude was estimated to be as large as eB = 10^''' 
Gauss ~ m^. Possible experimental signatures of the effects of this magnetic fleld have 
been suggested previously in conjuction with triangle anomaly; chiral magnetic/ separation 
effects IHl nil HI] and chiral magnetic waves [Ml [15] , leading to charge dependent two- 
particle correlations [16] and charge-dependent elliptic flows of pions [ITJ [18]. See also 
Ref.[T9]. These predictions are indeed observed experimentally at STAR/PHENIX 
[20| l2n [221 [23] and LHC [21] , and more experimental results are on the way. 

The idea that the early magnetic fleld may affect photon emissions leading to az- 
imuthal asymmetry has been previously discussed in Refs. [251 [SHI [271 12H1 [29] , mainly in the 
weak coupling regime [251 ES] or including conformal anomaly [27] . In our present study, 
we look at the effect of magnetic fleld on the photon creation, especially its azimuthal 
asymmetry that leads to observable mock-up of elliptic flow, in a strongly coupled flnite 
temperature plasma in the framework of gauge/gravity correspondence. See Refs. [30] [31] 
for earlier holographic computations for the photon emission rate without magnetic fleld, 
and Refs. [32] [33] for the photon emission with an anisotropy introduced by an external 
axionic perturbation [31]. Ref. [35] [36] computed photon emission rates with magnetic 
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field, but the detailed elliptic flow has not been computed. As the photon emissions are 
sensitive to microscopic details of the theory, we choose to work with the model by Sakai 
and Sugimoto [37j which is dual to the fleld theory which should be closest to the real 
QCD in quenched approximation. 

Our results indicate that the azimuthal dependence of the photon emission rates is 
more complicated than the simple elliptic flow pattern, which necessitates looking at 
higher moments of the azimuthal dependence, such as the quadrupole flow v^. We predict 
an interesting pattern of their ratios, j^^, as a function of frequency, which might be 
relevant experimentally. 

We also compute polarization asymmetry of the photons emitted in the presence of 
the magnetic fleld. After deflning an observable, "in/out plane polarization asymmetry", 
with respect to the reaction plane, we compute it as a function of frequency. Al- 
though it seems challenging to measure polarizations of photons in current experiments, 
it might be an interesting direction to reflne the experimental photon measurements. 

In both computations, we observe that triangle anomaly of the chiral symmetry plays 
an indispensible role; in gauge/gravity duality it manefests as a 5-dimensional Chern- 
Simons term which does affect the equations of motion in the presence of magnetic fleld. 
Essentially the effects can be summarized as coming from the new transport phenomenon 
of chiral magnetic waves of chiral charges. We explicitly check that results with and 
without Chern-Simons term differ substantially to each other, which indicates that any 
computation without properly including chiral anomaly might be questioned. 

In our analysis, we assume a static plasma and a constant magnetic fleld. In compar- 
ison to experiments, one has to integrate our results over the realistic time history of the 
heavy-ion collisions. This procedure is based on an additional assumption of adiabaticity. 
We leave further reflnement of our study including time-evolution of the plasma to the 
future. 

2 Photon emission rates from equilibrium plasma 

In this section, we summarize photon emission rate formula and useful relations between 
Green's functions for the case of translationally invariant thermal plasma. Note that we 
will not assume rotational isotropy that is broken by the presence of a magnetic fleld. 
The differential photon emission rate summed over all polarization states from a plasma 
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in thermal equilibrium is given by 



(2vr)3 2|fc 



(2.1) 



fcO=|A:| 



where is the Wightman function defined by 



The emission rate for a particular polarization state specified by is 



dr^ 



d^k 

(27r)3 2|^| 



(2.2) 



(2.3) 



fcO = |fc| 



Our metric convention is rj = (— , +, +, +). We would like to relate with the retarded 
Green's function G^ defined by 



G^Ak) ^ (-0 / d'xe-^'^9ix') J^^^(x), 7^(0) 



(2.4) 

since G^ is naturally computed in holography with incoming boundary condition on the 
black-brane horizon as we will see in the next section. It is convenient to define 



G>.(fc) 



d^xe-*'="(jf^^(a;)jf^*^(0)) , 



and using Lehmann representation one can easily derive the relation 

GUk) = -e^'''G<Ak) . 
On the other hand, from definitions of G^'^ one has 



(2.5) 



(2.6) 



G^Ak) 



d xe~ 



ikx 



^."""'(O)]) = -GUk) - G^k) = {e''' ' l) G<,{k) 



(2.7) 

Looking at the two definitions ( 12. 4p and (12. 7p . one has the following relation between G^ 
and G^, 



G^k) 



1 r+°° ~ 1 
27r i-oo k^-k^ - is 



-J-r d~k' ^ G^Ck^ k) - -GZ{k) 



{21 



where V denotes Cauchy principal integral. The relations ( 12. 7p and (12. 8p are the basic 
starting point of our discussion. 



From translational invariance and hermiticity of Jf'*^; easily derives 



= = , (2.9) 

so that as well as and are hermitian matrices with (/iz/) indices, which means 
importantly that 

ri^^Glik) , e^e^*Gl{k) , (2.10) 
are real. Contracting (12. 8p with t]^^ or e^e^* therefore gives one that 



Note that 



-2Im 



V^'^Glik)] , e^e^*GZ{k) = -2 Im {e^e-'^G^M 



Im 



(2.11) 



(2.12) 



in general, especially for circular polarization states where is complex. Note also that 



G;:,(A;)^-2Im G^^{k) 



(2.13) 



in general, as G is hermitian but may not be real. Only in the case of isotropic plasma, 
G^^ik) (and G^'^) is diagonal and hence real, then the equality in f l2.13p holds. 



From (12. 8 p and (12. lip , one finally has 



d^k 



{2Tif2\k\ e^fc" -1 



Im 



fcO=|A:| 



and 



dV^ {e^') 



d^k -2 



Im 



e'^e'^*G^Jk) 



kO=\k\ 



(2.14) 



(2.15) 



which we use for our numerical computation of photon emission rates. 



3 Computational set-up in holographic QCD 

For our non-perturbative computation of photon emission rates in the presence of magnetic 
field at strong coupling, we work in the Sakai-Sugimoto model; a holographic model which 
is closest to the large A^^^ QCD with massless chiral quarks in quenched approximation. 
The chiralities of massless quarks and the aspects of triangle anomaly manifest themselves 
in the model through 5-dimensional Chern-Simons terms in the action, which will play 
important roles in our results later. The model was constructed originally via embedding 
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D8/D8 branes inside the holographic background generated by Nc number of DA branes. 
Roughly speaking, Nc DA branes describe Yang-Mills gluonic dynamics of large Nc QCD, 
whereas D8/D8 branes describe massless chiral quarks of left-handed/right-handed chi- 
ralities respectively. Since electromagnetism couples to the vector-like global symmetry 
of these chiral quarks, we are interested only in the dynamics of the 8-branes in quenched 
approximation for our photon emission rate computation. The original model sits in a 
10-dimensional geometry which is a warped product of 5-dimensional holographic space 
and an extra internal space of S'^ x S^, but for our purposes in this paper the internal 
space doesn't play a meaningful role, and we will integrate out our action over it to have 
a reduction to the 5-dimensional holographic space from the start. 

The resulting D8/D8 brane action is 

Sds/d-8 = -Ci?t J d^xdUU'^^-det{glj, + 27rP^F) 

^ ^2 j d'xdU e''''''^''AMF^pFQn , (3.16) 



where e-tensor is numerical and 

N? 



C= TT-riK , R'^ngsNJl . (3.17) 



The sign of the Chern-Simons term in the second line encodes chirality that the 8-branes 
describe; the D8 brane with the upper sign describes the left-handed U{Nf)l chiral 
dynamics, and vice versa for D8 and the right-handed chiral symmetry U{Nf)r- In the 
deconfined phase of the theory that we are interested in here, the 5-dimensional geometry 
in Eddington-Finkelstein coordinate is given as 



dsl^^(^y (-f{U)df + ^{dx^y]+2dUdt , /([/) = l-(^y , (3.18) 



i=l 

which has a black-hole horizon at U — Ut with 



in terms of temperature T. Since the final gauge theory observables are independent of 
Is, one can conveniently choose 27rZ^ = 1 which will be assumed throughout this paper. 
Sakai and Sugimoto determined the values of parameters in the model by fitting to the 
pion decay constant and the p meson mass with Nc — 3. In units of GeV, this gives us 

C = 0.0211 , = 1.44 , (3.20) 



6 



which defines the model without further free parameters. 

In this background, the DS and 1^8 branes do not meet each other, and their dynamics 
are independent of each other in leading A^^ approximation. The relation between these 
chiral fields on the 8-branes and the more conventional vector/axial fields is given by 

Av,a = 1{±Al + Ar) , (3.21) 

where Al^r mean the gauge fields on D8 and DS branes respectively. Note that in terms 
of currents, this is equivalent to 

Jv,a^±Jl + Jr . (3.22) 

Electromagnetism couples to the vector symmetry with the coupling strength e, which 
means practically for us that we replace 

Av bAem , (3.23) 

where Aem is the electromagnetic gauge potential. In our situation of having constant 
magnetic field along say direction, this implies that we have 

F^^ = eB , Aa = , (3.24) 

which is equivalent to having the same background magnetic fields on both DS and DS 
branes, 

F^^^F^^^eB . (3.25) 

It is easy to check that these constant background magnetic fields satisfy the equations 
of motion of 8-branes. 

Working in the Eddington-Finkelstein coordinate which is useful for computing re- 
tarded Green's functions 

dsl^^i^" i^-f{U)de + j^{dx'f^+2dUdt , /([/) = l-(^)' , (3.26) 

and using the identity 

^det(l + M) = 1 + ^trM + ^ (trM)^ - ^tr (m^) + O (m^) , (3.27) 
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the quadratic expansion of the action in the presence of the background magnetic field 
Fi2^ = eB becomes 

Ci2) = \ (a{U) {F,uf - B{U) {F,uf - C{U) {{F,uf + {F^uf) - D{U) {{F,,f + [F^^] 
- E{U) (F12)' + 2F{U)Ft3F3u + 2G{U) {FnF^u + Ft2F2u) j 

T {AtF^u - A^Ftu + AuFt^) , (3.28) 

where the coefficient functions are given by 



A{U) = CU^U^ + {eBYR^ , B{U) = CU f{U)^U^ + [eBfR^ 



EiU) = C(^)' — J , FiU) = CU (^) ' ^U^ + {eBfR^ , 

G{U) = Ci-] , . (3.29) 



du{B{U)F3u)-D{U){diF,3 + d2F23) + F{U)dtF3u-du{F{U)Ft3)T^^Ftu = 



The equations of motion derived from the above are 

A{U)dtFtu - B{U)d^F^u - C{U) {d,F,u + d^F^u) + F{U)dA 

N eB 

+GiU)idiFa + d2Ft2)±^;^Fr, = 

du {A{U)Ftu) + F{U)dsFsu + G{U) {d^F^u + ^2^21/) T ^^P^u = 

N^eB 

du{C{U)F^u) + D{U)d;,F^:, + E{U)d2F^2 + G{U)dtF^u-du{G{U)Ft^) = , 
du{G{U)F2u) + D{U)dA-E{U)diF^2 + G{U)dtF2u-du{G{U)Ft2) = . 

(3.30) 

Note that triangle anomaly represented by the 5D Chern-Simons term contributes to 
the equations of motion. Physically, this is due to an interplay between vector and 
axial symmetries via chiral magnetic waves in the presence of the magnetic field. This 
implies that the resulting photon emission rate encodes non-trivial effects from the triangle 
anomaly. 

We work in the gauge Au = 0, and solve the equations of motion (I3.30p for A^{U) 
in the frequency-momentum space assuming the factor e^''^ = Because there 



8 



is a residual 5*0(2) rotation symmetry on the transverse coordinates {x^,x^), we can put 
= without loss of generality. As we are interested in the elliptic flow, we write 

A;^ = A;°sin^ , k^^k° cos 6 , (3.31) 

satisfying the on-shell condition = \k\ = co and showing the azimuthal angle 9 from 
the direction of the magnetic field explicitly. The equations of motion then become 

A{U)duAt + cos 9B{U)duA^ + sin 9C{U)duAi - ioo cos 9F{U) {A3 + cos 9 At) 

N eB 

-iujsm9G(U)(Ai + sm9At)^ ^^(As + cos9At) = , 

N eB 

du{A{U)duAt)+iujcos9F{U)duAs + iujsm9G{U)duAiT^^duA3 = , 

du {B{U)duA3) - uj^ sin 9D{U) {sm9A3 - cos 9 A^) - iuF{U)du A3 

N eB 

-iujdu{F{U){A3 + cos9At))^^-^duAt = , 

du (C(U)duAi) + u'^ cos 9 D(U) (sinMg - cos Mi) 

-iujG{U)duAi-iujdu{G{U){Ai + sm9At)) = , 

du {C{U)duA2) - cos^ 9D{U)A2 - uj^ sin^ 9E{U)A2 

-iujG{U)duA2-iujdu{G{U)A2) = . 

(3.32) 

As we arc interested in the retarded Green's function of the global symmetries that the 
above 5-dimensional gauge fields describe, we need to establish a well-defined procedure 
of extracting information on the corresponding currents in the field theory side from 
the 5-dimensional profile of the solutions of the above equations. We follow standard 
steps of gauge/gravity duality dictionary. We first discuss how to extract expectation 
values of the symmetry currents in the field theory side from the 5-dimensional solutions. 
To construct the boundary current expectation value (J^) from the solution Af^{U) of 
the above equations, one has to study the near boundary asymptotics at f/ — ?> 00 and 
perform a careful holographic renormalization. The necessity of a careful holographic 
renormalization comes from the fact that we will turn on space-time varying external 
sources to find retarded Green's functions at finite frequency-momentum. In this situation, 
it generally happens that the renormalized boundary currents (J^) are given not only 
by the coefficients of the subleading terms in the asymptotic expansion, but also get 
some additional contributions from the external sources. Only after combining these two 
contributions, the resulting currents (J^) satisfy the correct conservation Ward identity. 
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The near U ^ oo asymptotics of are found to be 

U 2 U U 2 U 2 

where we assume the dependence e*'^'^ = e^*^°*+*'^'^ imphcitly. We will consider a general 
/c^ for a while in our discussion of holographic renormalization. A^^^ is the external source 
and the subleading terms A^^\ A^^\ and A^^^ are completely fixed by the external source 
A^^^ in the following way, 

4') = o , 4^) = , 

Af) = -2ii?i (A;Mf ) + Mf) = 2i?iFf , i = l,2,3 , 
aI'^ = -2R'{k'k^Af + k^k,A?^) = -2R%F,f , 
Af^ ^-2R^(k^kjAP -kik^Af) = -2R%F^^ . (3.34) 

The fact that the coefficient of the logarithmic term, ^4^^^ vanishes is special in this model, 
which would not be the case for a 5D gauge theory in asymptotic AdS^ geometry. The 
absence of logarithmic term translates to the absence of contribution to the conformal 
anomaly from the corresponding global symmetry in the 4D gauge theory side. This can 
be understood because the 4D theory that this holographic model supposedly describes 
is conformally non-invariant even in ultra-violet regime. The piece A^ is not determined 
by the external source, and it encodes a dynamical freedom of the expectation value {J/j). 
It should be fixed by appropriate infrared boundary conditions, which in our case will 
be the incoming boundary condition on the black hole horizon. In Eddington-Finkelstein 
coordinate we are working, this simply requires regularity of A^{U) at the horizon U — 
Ut- In conjunction with the ultraviolet boundary condition A^^\ these two boundary 
conditions determine the solution Afj_{U) uniquely. 

Although {Jfj} should contain A^ representing dynamical degrees of freedom, the full 
expression for (J^) involves additional contributions from (derivatives) of the external 
source A^^^ as we mentioned above. To find these additional contributions, one normally 
goes through steps of holographic renormalization by regularizing and subtracting diver- 
gences as we describe in the following. For this purpose, it is more convenient to go to 
the diagonal metric frame. 



i=l 
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by the coordinate transformation, 



oo 



Near f/ — t- oo boundary, the temperature/magnetic field do not matter and there is 4D 
Lorentz symmetry one can use in discussing asymptotics and expectation values in this 
frame. The near boundary behavior of A* in the AIj = gauge in this frame is found to 
be (we put * to mean the gauge fields in this diagonal metric frame) 

A; = A(^^~2R'd.F('>^ + ^ + ... . (3.37) 

u (7 2 

Note the absence of U^^ term as well as the logarithmic term. Upon inserting the above 
to the 5D holographic actioij*] 

= - J / d'^dU [2UiF^uF^'' + ^J^.F^"^ , (3.38) 

there appears a divergence, 

SId ~ -^R'ui I , (3.39) 

where we regularize the divergence at f/ = Uoo- The counterterm at f/ = Uao is hence 
chosen to be 



= +-R^ui J d^xF^.F^'^ 



(3.40) 

to cancel the divergence minimally. Then the renormalized current expectation value from 

gren ^ greg ^ get 

W = lim (-CukduA^{Uoo) + 2CR'ukd,F^''{U^)) = lcA* . (3.41) 

The result is simple without further contributions from the source; this is special in this 
model, which is related to the absence of logarithmic term in the expansion f l3.37p . 

The 5D Chern-Simons term does not introduce further divergences, so the minimal 
counterterm (I3.40p remains unchanged. It only modifies the regularized action S^*^^ by 
(recall our A^ = gauge) 

55^^^ = I d^xdUe'^'^-'A^F^u d^xdU e^"''^ A^F^^Fpu , (3.42) 



*One can neglect effects from temperature and magnetic field in near [/ — )■ oo boundary as they 
are sufficiently subleading. The 5D Chern-Simons term is also irrelevant in divergences and minimal 
counterterms and we will consider it only at the end in the final results. 
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with e*^^"^ = 1, and is the gauge field responsible for the background magnetic field, 
Afg = =F^x^'^. The correction to (J^) from this is thus 

W = T^e,,,.A^'> ± . (3.43) 

When we consider both D8 and D8 branes together in discussing vector/axial symme- 
tries, one can still add a further finite piece in the counterterm S"*^*, called the Bardeen 
counterterm, to ensure the conservation of vector symmetry, and this will modify (J^) 
additionally [38] . 

Although one may deal with them carefully to get the final results correctly, for our 
problem of computing vector- vector Green's functions it turns out that these modifications 
from f l3.43p and the Bardeen counterterm simply drop in the final results. Recall that the 
vector current is a simple sum of left- and right-handed currents, so that the modification 
f l3.43p gives rise to 

where A^^ = | (— + ^^^) is the source for the axial current Ja = ~Jl + Jr- The 
modification from the Bardeen counterterm is also easily found to be 



X/ T \ _ NceB a(^> , aBu p(0)o^P foA^\ 

6Ti LZTX^ 



Both (13.441) and (13.451) give the response of ( Jy) proportional to the external axial source 
Aa only, so that they don't contribute to the vector-vector Green's functions. 

Based on this observation, we will use the formula (I3.4ip for the current expectation 
values for simplicity, without worrying about further corrections from the 5D Chern- 
Simons term. We stress that the presence of the 5D Chern-Simons term does affect 
our results via bulk equations of motion (13.301) : it affects the dynamical coefficient 
and hence the current expectation values through equations of motion. This can easily 
be seen by adding and subtracting (I3.30p for upper/lower signs of Chern-Simons term 
(corresponding to D8/D8 branes for left /right-handed chiralities) to get 



A{U)d,F^/ - B{U)d,F^^ - CiU) {dXif + d2Fi,^) + F(U)d,F, 



,V,A 
t3 



-fG(f/)(5,Fr + 5.Fr)-^F.r = 



An' 

" Ay 



du {A{U)Fru^) + F{U)d,F^^ + G{U) {d^F^^ + d,F^^) + ^F^u 



= 
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Figure 1: Contribution of triangle anomaly to the vector- vector Green's functions in the 
presence of magnetic field. 



du {b{U)F^^) - D{U) [d,F^^ + d,F:^i^) + F{U)d^F, 

-du [F{U)Fr^ 



V,A 
3U 



+ 47r2 



du {C{U)F^^) + D{U)d,F^^ + E{U)d2F'^^ + G{U)dtF^{f - du {G{U)F, 
du {C{U)F^lf^) + D{U)d,F^f - E{U)d,F^'' + G{U)d,F^^'' - du {G{U)F, 



V,A 
tl 



V,A 
t2 



, 

, 
, 
(3.46) 



where A^'^ = | {^A^ + Aji). It is clear that axial components are necesssarily excited 
in computing vector-vector Green's functions due to Chern-Simons term. A diagramatic 
representation is given in Figure [TJ Therefore, the effects from triangle anomaly to our 
photon emission observables are totally dynamical, and are not sensitive to the issues of 
additional boundary contributions. In retrospect this makes good sense because at the 
end we are dealing with only vector-like observables of electromagnetism which should be 
unambiguously defined. 

What remains to find (J^) in the original Eddington-Finkelstein coordinate is to per- 
form the coordinate transformation of ( 13.36^ on A^{U), so that one can express A* in 
terms of and A^^\ One has to be careful about our gauge choices Au = and = 0; 



13 



starting from A^{U) and performing the coordinate transformation fl3.36p . the resulting 
is not zero and one needs to do a further gauge transformation to remove AIj. Note 
that dt = df , but du ^ du* ■ The computation is straightforward and the result is 



Al = A + -Rld,d,F,f 



,9 



A* = A + AR^^[d,d,F}^^ + ^d^F}p-^d,d,F,f) , (3.47) 

from which one obtains the current expectation values via (J^) = as in (13.41 p . One 

can check that the first equation in (I3.30p gives one the correct chiral Ward identity (up 
to the additional boundary contributions and the Bardeen counterterm discussed above) 

MJn=±^FS' , (3.48) 

in the presence of external chiral gauge fields eB and A^^\ 

Once we know how to extract the expectation values from the solutions as above, 
we can easily compute the retarded Green's functions in the following way. The leading 
component A^-^^ in the near boundary expansion ( I3.33P is interpreted in the QCD field 
theory side as an external gauge potential coupling to the (chiral) current r- The 
chirality depends on the sign of the Chern-Simons term, or equivalently on whether we are 
looking at D8 or D8 brane. By demanding incoming boundary condition at the horizon, 
which is simply a regularity at the horizon in our Eddington-Finkelstein coordinate, the 
solution is uniquely determined by this external source ^4^°^; it is clear that the solution 
linearly depends on A^-;'\ One then obtains the expectation value of chiral current from 
the solution by using (I3.4ip and (I3.47p . As the result is linear in A^^\ one writes 

W = -G^'{k)A^^^ , (3.49) 

which gives one the retarded Green's function G^^^ by the definition of Kubo's formulation 
of real-time response functions. For the upper sign of Chern-Simons term in (I3.30p ( that 
is from D8 brane), one obtains the Green's function of left-handed chiral current G^^, and 
vice versa for the lower sign (that is D8 brane) and the right-handed G^^. The desired 
vector-vector Green's function for the photon emission rate is then 

Gyv = + G^j^ . (3.50) 

Because the two chiral Green's functions are simply related to each other by 

GLLieB) = GRRi-eB) , (3.51) 
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or equivalently 

GLL{0) = GRR{n-e) , (3.52) 

in terms of azimuthal angle 6 with respect to the magnetic field, it is enough to compute 
Gll only. 

We numerically solve the equations of motion f l3.30p . or more precisely (13.321) . to 
compute the retarded Green's functions from the relation (I3.49p . Instead of solving the 
equations with a given A^^^ from the ultraviolet (UV) boundary U — cxd, it is numerically 
much easier to solve them from the horizon U = Ut toward the UV boundary, starting 
with a regularity boundary condition at the horizon. Inspecting the equations of motion 
(I3.32P near the horizon U = Ut-, the regularity condition uniquely fixes the derivatives 
duA^iUr) and dfjA^iUT) in terms of the value at the horizon A^{Ut), which allows one 
to start solving the equations (I3.32p from the horizon. Practically, we use djjA^{UT) 
and dijAfj^iUx) to proceed a small step to U = Ut + ^ with a small number e = 0.01, 
and start our numerical solving the equations ( 13.32^ from that point until a UV cutoff 
U = Umax- Given this solution, one obtains A^^^ and A^ by comparing with the near 
boundary expansion (I3.33P and (I3.34p at the position U = Umax- More precisely, we solve 
the following linear system of equations for {A^^\ A), 

/1(2) A /1(4) 

A (U ) = v4(°) + ^ + ^ + ^ + ^' 

^fiy-^max) ' I ' tt ' S. ' TT2 ' 

Tn 'J max Tn 'J ^ 

Umax Umax 



max 



1 /l(l) /1(2) o A A{A) 

{duA,){Umax) = -i^- 77^ -^4^ - 277^ ' (3-53) 



/max 



where the left-hand sides are given by the numerical solution, and A^^'^'^^ are linearly 
given by y4^°) as in (I3.34p . Note that we have included one more term, A'^'^\ in the 
above expansion on the right-hand side than in (I3.33P for a better numerical precision of 
extracting {A'^^\ A) from the above. The expressions for A^^^ can be easily found and we 
skip their explicit expressions. 

From (A*^°\ A), one then constructs {A^^\ (J)) using (I3.4ip and (I3.47p . Since it is clear 
that the results are linearly dependent on the horizon data A^{Ut) that we start with, 
one writes 

Af =S^''A,{Ut) , {J,)=n;A,{UT) , (3.54) 

with two matrices S and TZ- One easily computes these two matrices by performing the 
above described numerical procedures for each four unit vectors of A^j_{Ut)- Once they 
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are found, one relates {A^^\ (J)) directly by 

W = 7^/ (^5-0/^^°^ = ■ ^5-0/ 4°^ - -G^ ^^1°^ , (3.55) 

so that the retarded Green's function is finally computed as = —TZ ■ in a matrix 
form. In the total photon emission rate formula f l2.14p . one needs the trace of the Green's 
function and 

V'^Gl = Gy = -tr (7^ ■ S-') , (3.56) 

which is particularly simple in this matrix form. However, the full matrix Green's function 
contains much more information than the simple trace, such as polarization states of 
emitted photons. We will discuss polarization asymmetry of the photons in section [5l 
which we propose to be an important new experimental signature of triangle anomaly in 
heavy-ion experiments. 

From explicit expressions for A^^''^'^\ one notes that the expansion parameter in the 
near boundary series (I3.53P is 

X = , (3.57) 

U max 

SO that the larger the energy u is, the bigger Umax should be for a good numerical precision 
of A^^^ and (and hence of the retarded Green's functions) from solving the equations 
f l3.53p . From f l3.53p . the error on {A^^\ A) from neglecting higher order terms on the right- 
hand side is of order x^. In our numerical analysis, we have controlled to be much less 
than 0.01, so that our numerical results should be reliable up to better than 1%. Also, 
the structure of equations of motion f l3.30p is such that the derivative of the first equation 
in (I3.32P is equivalent to the remaining four equations. Therefore, once the first equation 
is satisfied at one point in U, say at f/ = Ut, it should remain satisfied for all U if one 
solves the other equations correctly. In other words, it is a consistent contraint equation 
of the system. It serves in fact useful in our numerical analysis by providing a test of the 
numerical precision, and we have checked that our solutions satisfy it with much better 
precision than 1%. Another independent check is the chiral Ward identity fl3.48p which 
also holds to a good precision. This gives us the conservation Ward identities of our final 
vector-vector Green's function, 

FG^K -^(k) = ''{k)K = , (3.58) 

which are checked to be true with a very good precision. 
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Figure 2: Azimuthal angle dependence of photon emission rates for three different en- 
ergies u = 0.1,0.5, 1 GeV with temperature T = 0.2 GeV and magnetic field eB = 0.4 
GeV^. The dashed curves are the results when the Chern-Simons term is turned off. 

4 Elliptic (^2) cind quadrupole (^4) flows of photons 

With the computations described in the previous section, we present our results for the 
azimuthal dependence of the photon emission rates in the presence of the magnetic field. 
We introduce the mode expansion in the azimuthal angle 6 with respect to the magnetic 
field direction as 

^(^) = ro(l-2t;2cos(2e) + 2t;4cos(4^) + ---) , (4.59) 

with the elliptic flow V2 and the quadrupole flow v^. The negative sign in front of V2 in 
the above definition is due to the relation 

^=^-0 , (4.60) 

between 6 and the more conventional angle from the reaction plane (note that magnetic 
field is perpendicular to the reaction plane). In Figure [21 we show some exemplar plots 
of azimuthal dependence of photon emission rates ^^(^) for different energies. We take 
T = 0.2 GeV, eB = 0.4 GeV^ for an illustrative purpose (we will also consider more 
realistic values of the magnetic field later). For useful comparison, we also include results 
obtained after dropping the Chern-Simons term (the dashed curves), which shows the 
importance of triangle (chiral) anomaly in the results, especially for low energy u < 1 
GeV. As is clear from the plots, the angular dependence is drastically affected by the 
triangle anomaly for low energy regime (the plots in the left), and the mode expansion 
becomes more complicated than simply being characterized by an elliptic flow. Especially, 
the size of the quadrupole moment is comparable in this regime. We expect that these 
modifications of the azimuthal dependence are due to the existence of the chiral magnetic 
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Figure 3: Elliptic flow V2 versus photon energy u for T = 0.2 GeV and eB = 0.4 GeV^. 
The right plot is the result without Chern-Simons term (triangle anomaly). 



wave modes that affect the retarded Greens functions via its pole structure 

1 1 

+ 



(CO 



f^/ccosl 



[u + v^k cos I 



(4.61) 



where is the velocity of the chiral magnetic wave and the two contributions are from 
left- and right-handed chiral magnetic waves respectively. We plot f 2 (in Figure [3]) and 
the ratio f 4/(^2)^ (in Figure Hj) to highlight this effect, for example the violation of the 



usual scaling 



In the elliptic flow V2, we notice that the elliptic flow can even 



be negative for the low energy u < 0.3 GeV due to the effect from triangle anomaly. 
Current experiments give data only for u > 1 GeV, so that the experimental relevance 
of this observation is not high at the moment, but it may become important in the 
future. In Figure O we plot the elliptic flow V2 for a realistic value of the magnetic 
field eB = 4m^ ^ 0.08 GeV^. The general trends are the same as we describe in the 
above, although we observe that the overall magnitude of the elliptic flow is rather small, 
V2 ~ 10^^. We should take this as a result from the strong coupling computation via 
gauge/gravity correspondence. The particular reason why the strong coupling gives such 
a small azimuthal imbalance of the photon emission even with a large strength of the 
magnetic field is not completely clear to us, and might deserve a further study. 



5 Polarization asymmetry of photons 

Our computational framework in gauge/gravity correspondence is able to describe the 
photon emissions with specific polarization states of the photons. With our previous 
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Figure 4: The ratio v^/ {V2Y versus photon energy a; for T = 0.2 GeV and eB = 0.4 GeV^. 
The right plot is the result without Chern-Simons term (triangle anomaly). 



T=0.2 GeV, eB=0.08 GeV^ 



0.2 0.4 • 0.6 0.8 



Figure 5: Elliptic flow V2 versus photon energy u for T = 0.2 GeV and eB = 4m^ = 0.0^ 
GeVl 



choice of photon momentum 



= Lo (1, sin6', 0, cos( 



we deflne the in-plane polarization with 

£^^ = (0,0,1,0) , 
and the out-plane polarization with 

^OUT — (0; "^OS ^; 0; ~ ^in ( 



(5.62) 



(5.63) 



(5.64) 



See Figure M for a schematic explanation of our deflnition of the polarization states. We 
then deflne a new observable, the "in/out-plane polarization asymmetry" A^^'^ by 



4I/0 _ J3 



^(cin) + ^(cout) 
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(5.65) 




1 

X 



Figure 6: Definition of in- and out-plane polarizations. 



T=0.2 GeV. eB=0.08 GeV-, iO=Q.i GeV t=0.2 GeV. eB=0.08 GeV-, w=0.5 GeV T=0.2 GeV. eB=0.08 GeV', w=l GeV 




Figure 7: In/out-plane polarization asymmetry A^^'^ for three different energies oj = 
0.1,0.5,1 GeV with temperature T = 0.2 GeV and mag netic field eB = Ami = 0.08 
GeV^. The dashed curves are the results without Chern-Simons term (triangle anomaly). 

In Figure [71 we plot the azimuthal dependence of the in/out-plane polarization asymmetry 
j^i/o fQj, different energies. As before, we observe that the Chern-Simons term affects the 
results in a fundamental way for low energy regime u < 1 GeV, and we expect that the 
chiral magnetic wave pole to be responsible for this modification. The general trend is 
that A^^'^ is positive for low energy regime to < 1 GeV, whereas it becomes negative for 
higher energies u > 1 GeV. Although current experiments seem to find it difficult to 
produce data for the photon polarization, our observation may become relevant in more 
refined experimental measurements in the future. 
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